Control of Stationary Behavior in Probabilistic Boolean Networks by Means of Structural Intervention

Probabilistic Boolean Networks (PBNs) were recently introduced as models of gene regulatory networks. The dynamical behavior of PBNs, which are probabilistic generalizations of Boolean networks, can be studied using Markov chain theory. In particular, the steady-state or long-run behavior of PBNs may reflect the phenotype or functional state of the cell. Approaches to alter the steady-state behavior in a specific prescribed manner, in cases of aberrant cellular states, such as tumorigenesis, would be highly beneficial. This paper develops a methodology for altering the steady-state probabilities of certain states or sets of states with minimal modifications to the underlying rule-based structure. This approach is framed as an optimization problem that we propose to solve using genetic algorithms, which are well suited for capturing the underlying structure of PBNs and are able to locate the optimal solution in a highly efficient manner. Several computer simulation experiments support the proposed methodology.